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ABSTRACT 

Aims. We aim at characterizing a sample of nine new hard X-ray selected Cataclysmic Variable (CVs), to unambiguously identify 
them as magnetic systems of the Intermediate Polar (IP) type. 

Methods. We performed detailed timing and spectral analysis by using X-ray, and simultaneous UV and optical data collected by 
XMM-Newton, complemented with hard X-ray data provided by INTEGRAL and Swift. The pulse arrival time were used to estimate 
the orbital periods. The broad band X-ray spectra were fitted using composite models consisting of different absorbing columns and 
emission components. 

Results. Strong X-ray pulses at the White Dwarf (WD) spin period are detected and found to decrease with energy. Most sources are 
spin-dominated systems in the X-rays, though four are beat dominated at optical wavelengths. We estimated the orbital period in all 
system (except for IGR J16500-3307), providing the first estimate for IGR 108390-4833, IGRI18308-1232, and IGRJl 8 173-2509. 
All X-ray spectra are multi-temperature. V2069 Cyg and RX J0636+3535 posses a soft X-ray optically thick component at kT~80 eV. 
An intense K„ Fe line at 6.4 keV is detected in all sources. An absorption edge at 0.76 keV from OVII is detected in IGR J08390-4833. 
The WD masses and lower limits to the accretion rates are also estimated. 

Conclusions. We found all sources to be IPs. IGR J08390-4833, V2069 Cyg, and IGR J16500-3307 are pure disc accretors, while 
IGR 118308-1232, IGR J1509-6649, IGR J17195-4100, and RX 10636+3535 display a disc-overflow accretion mode. Afl sources 
show a temperature gradient in the post-shock regions and a highly absorbed emission from material located in the pre-shock flow 
which is also responsible for the X-ray pulsations. Reflection at the WD surface is likely the origin of the fluorescent iron line. There 
is an increasing evidence for the presence of a warm absorber in IPs, a feature that needs future exploration. The addition of two 
systems to the subgroup of soft X-ray IPs confirms a relatively large (~ 30%) incidence. 

Key words. Stars: binaries: close - Stars: individual: IGR 108390-4833, IGR 118308-1232, IGR J16500-3307, IGR J18173-2509, 
IGR 117195-4100, V2069 Cyg (also known as RX I21237-I-4218), RX 10636+3535 (also known as V647 Aur), IGR 115094-6649, 
and XSS J0056+4548 (also known as V515 And) - X-rays: binaries - Accretion, accretion discs - stars:novae, cataclysmic variables 



1. Introduction 

Magnetic Cataclysmic Variables (CVs) constitute a subgroup of 
the CV class, harboring accreting white dwarfs (WD) with mag- 
netic field strengths B ^ 10^ G. These systems are further sub- 
divided in two groups, depending on the WD magnetic field in- 
tensity and degree of synchronism (Prot=M/Poib=n)- Those called 



* Based on observations obtained with XMM-Newton and 
INTEGRAL, ESA science missions with instruments and contri- 
butions directly funded by ESA Member States; and Swift, a NASA 
science mission with Italian participation. 



Polars are synchronous mCVs. They show signatures of strong 
magnetic fields (B ~ 10 - 230 MG) through the presence of 
a conspicuous polarization at optical and near-IR wavelengths. 
The so-called Intermediate Polars (IPs) possess instead asyn- 
chronously rotating WDs (P„ << Pn) and except for a few cases 
(Piirola et al. 1993 ; Potter etal. 1997; Katajainen et al. 20071 
iButters et al.ll2009l; iPotteret al.ll2012l) . they do not show polar- 
ization. These systems are believed to possess weakly magne- 
tized WDs (B <10 MG). The Polars populates the orbital period 
distribution at short periods mainly below the 2-3 h orbital pe- 
riod "gap" (.Warner 1995: Wheatlev 1995.) . On the other hand, 
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IPs are generally found above the gap, with only a handful of 
systems below it. 

The material lost from the Roche-lobe overflowing late- 
type secondary star (typically a Main Sequence or a sub-Giant) 
may flow toward the WD in different modes, depending on the 
magnetic field intensity: directly (stream), through a truncated 
disc or ring. The high field Polars are direct accretors, while 
IPs may accrete via a stream or a disc/ring depen ding on their 
magnetic m oment and spin-to-orbital period ratio (iNorton et al.l 
12004 1200 8h. A combinatio n of these, called disc-overflow, may 
also occur (Hellier 1995; No rton et al.ll 19971 ) and is linked to the 
mass accretion rate. Close to the WD, the magnetically con- 
fined flow of matter is accreted through a colum n in stream- 
fed s ystems and through an arc-shaped curtain (iRosen et al.l 
Il988 l) in disc/ring-fed systems. Approaching the WD surface, 
the velocities are supersonic and a stand-off" shock is formed. 
In the post-shock region (PSR) the gas at a temperature of 
kT~ 10-60 keV slow s and cools via thermal bremsstrahlung and 
cyclo tron radiations (lAizulll973l: IWu et al.l[T994t ICropper et aP 
Il999l) . The efficiency of the two cooling mechanisms depends 
on the WD ma gnetic field strength (Woelk & Beuermann 1996; 
iFischer & Beuermannil200lh . The post-shock emission is par- 
tially intercepted and thermalized by the WD surface, giving rise 
to an optically thick emission that emerges in the soft X-rays 
and EUV/UV regimes. An intense soft (kT ~ 20-50 eV) com- 
ponent is a characteristic signature of the X-ray spectra of th e 
Polars, mostly balancing the cyclotron flux (lBeuermannll2004l) . 
On the other hand IPs have typically stronger hard X-ray fluxes 
than the Polars and only recently a soft h ighly absorbed opti- 
cally thic k component has been detected (de M artino et al.l2004t 
lEvans & HelUe r 2007; Anzolin et al. 2008, 200^^ 

Though mCVs represent a relatively small (~ 20%) frac- 
tion of CVs, this number is rapidly increasing thanks to the re- 
cent hard X-ray surveys conducted by INTEGRAL and Swift 
above 20 keV. S o far 64 CVs are identified in the latest cat- 
alogue releases dBird et al.l 1201 Ot ICusumano et al.l 1201 Oh . 43 
of them are mCVs. The IPs represent ~ 80% of this hard 
X-ray sample (including both already known (Barlow et al.' 
2006|) and recently discovered members (Bonnet-Bidaud et al. 
2007; de Martino et al. 2008; Anzolin et al. 2009; Butters et al. 
2008.; , Bonnet-Bidaud et al...2009; .Pretorius. .2009; Scaringi et ah 



201 1 ), though one misidentification was found (de Marti no et al.l 



2Qm . Hence, while in the pre-INTEGRAL and Swift times IPs 
amounted to only ~20 members, this number has more than dou- 
bled nowadays. 

Hard X-ray mCVs have the potential to be important con- 
tributors to the X-ray source population at low luminosities 



(~ 10 - 10 erg s ). T hey were proposed to be the maj or con- 
stitue nt of ga lactic ridge (ISazonov et al. 2006; Revnivtsev et al.l 
20081 120091) an d galactic bulge (Revniv tsev et all l201lt 



Hong et al. f l20T2h X-ray emission. Therefore, they are also be 



lieved to have an important role in the X-ray luminosity function 
of other galaxies. 

In the framework of a program conducted with 
XMM-Newton, with the main goal to identify new mag- 
netic type CVs, we present here the analysis of simultane- 
ous X-ray, UV and optical data of a sample of nine CVs 
that were detected as hard X-ray sources by INTEGRAL, 
RXTE or SWIFT: IGR J08390-4833, IGR J18308-1232, 
IGR J18173-2509, IGR J17195-4100, V2069 Cyg (also 
known as RX J21237-H4218), V647 Aur (also known as 
RX J0636+3535), IGR J15094-6649, V515 And (also known as 
XSS J0056H-4548), and IGR J16500-3307 (hereafter IGR J0839, 
IGR J1830, IGR J1817, IGR J1719, V2069 Cyg, RX J0636, 



IGR J1509, XSS J0056, and IGR J1650 respectively. We com- 
plement the X-ray analysis with a high-energy coverage from 
INTEGRAL/IBIS and Swift/BAT publicly available data. Five 
of these sources (IGR J0839, IGR J1817, XSS J0056, IGR J1509 
and IGR J 17 19) have also Chandra, RXTE, or Swift/XRT 
coverage. Afl, except IGRJ0839, IGRJ1830 and IGRJ1817, 
were also observed in optical photometry and spectroscopy 
from ground-based telescopes, from which orbital periods 
were determined and optical pulses detected, thus allowing a 
comparison with our results. 

In Sect. 13.11 we present the X-ray, UV and optical timing 
analysis. In Sect. l3.2l we present the analysis of the X-ray broad- 
band spectra of each source. For the three brightest sources 
we also present the high resolution spectra provided by the 
XMM-Newton RGS instrument. In Sect. |4] we discuss the ori- 
gin of rotational pulses, the accretion mode and the emission 
properties of our sample. 



2. Observations and data analysis 

2.1. XMM-Newton observations 

The XMM-Newton Observatory includes three ~ 1500 cm^ 
X-ray telescopes with an EPIC instrument in each focus, 
a Reflecting Grating Spectrometer, RGS (den Herder et 
,2001.) and an Optical Monitor, OM (Mas on et al. 200 iT 
Two of the EPIC imaging spectrometers use MOS CCDs 
(iTurneret al.l 1200 ll) and one uses a PN CCD (IStriider et al.l 
120011) . XMM-Newton coUected data of all nine CVs. The main 
observation parameters, together with that of all other observa- 
tory instruments are reported in table |T] 

Data were processed with SAS version 10.0.0, using the 
updated calibration files (CCF) available in January 2011. All 
observations were performed with the EPIC-PN (PN hereafter) 
camera set in prime full window imaging mode (time resolu- 
tion=0.0734 s) with a thin filter applied. Standard data screening 
criteria were applied in the extraction of scientific products. We 
accumulated a one-dimensional image and fitted the ID photon 
distribution with a Gaussian. Then, we extracted the source pho- 
tons from a circular region of radius 40" centered at the Gaussian 
centroid. The background was obtained (within the same CCD 
where the source lies) from a circular region of the same size. 
Single and double pixel events with a zero quality flag were se- 
lected for the PN data. For the spectral analysis we cleaned all 
observations from solar flares by collecting CCD light curves 
above lOkeV and applying an intensity threshold. Spectra were 
then produced only for those parts of observations under the 
threshold limit. On the other hand, for the timing analysis we 
used the epiclccor task, producing a background subtracted 
light curve in the range 0.3-15 keV (with a bin time of 1 1 s), con- 
sequently removing contamination from solar flares. The source 
event arrival time of each observation, in the 0.3-15 keV energy 
range, were converted into barycentric dynamical times (BDT) 
by means of the SAS tool barycen. 

The PN spectra were rebinned before fitting, to have at least 
30 counts per bin. We report the analyses obtained with the PN 
data only (consistence with the results of EPIC-MOS cameras 
was always verified). Phase-resolved spectra were also extracted 
at the pulse maximum and minimum. All spectra were analyzed 
using the version of XS PEC (12.5. In). 

We also inspected the RGS spectra. We used the 1st order 
spectra and responses produced by the pipeline, which runs an 
automatic version of the RGS meta-task, rgsproc. While, in 
principle, the quality of data products can be improved some- 
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what by manually repeating the reduction, we found no obvious 
problems with the pipeline products. Also considering the lim- 
ited statistical quality of the RGS data for our sample, we believe 
our choice to use the standard pipeline products is well justified. 
The RGS specti-a of IGR J 1 7 1 9, IGR J 1 509, and XSS J0056 were 
of high enough statistical quality to merit further analysis; these 
are discussed below in sect. 13.2.31 The data for the other targets 
were of too low S/N to add useful information. We grouped the 
RGS spectra of the three sources such that each new bin had at 
least 25 counts. 

Simultaneous coverage to the X-rays was ensured by the 
OM operated in fast window mode using sequentially the B 
(3815-4910A) or the U (3000-3800 A) filter and UVM2 (2000- 
2800 A) or UVWl (2450-3200 A) filter Consequently, each tar- 
get was observed for half of the time of the EPIC exposure in 
the B or U filter and the other half in a UV filter Due to the 
OM fast window mode acquisition, the OM photometric data 
in each filter consist of sequential (four or five) segments of 
length ~ 1 300-2800 s each (according to the total EPIC expo- 
sures). Average net count rates and instrumental magnitudes are 
reported for those sources with usable data in table[T] where we 
also report whether no variability is detected. Background sub- 
tracted OM light curves were obtained using the standard SAS 
processing pipeline with a binning time of 10 s. Barycentric cor- 
rections were also applied to the OM light curves. 



3.1.1. X-ray variability 

We computed the power spectra in the 0.3-15 keV range using 
the PN data. In all of them, strong peaks are found and in many 
cases also harmonics (see figure[ri|2] and|3]l. The periods of the 
main peak s were determin ed by means of a phase-fitting tech- 
nique, see iDairOsso et al.l (f2003 ) for details on the technique. 
The derived values are reported in table |2] All uncertainties are 
hereafter at Icr confidence level (c.l.) if not otherwise specified. 
For IGR J 17 19 we did not detect significant periodic variability 
in the last third of the observation and therefore, the results refer 
to the first 2/3 only. This source also shows a variability on time 
scale of a few thousands of seconds occurring at about 1/3 of 
the total observation time. This variability is likely responsible 
for the low frequency peaks seen in the power spectrum. We are 
however unable to characterize it with the present data. 



2.2. The INTEGRAL observations 

The INTEGR AL IBIS/ISGRI instmment (lUbertini et al.ll2003t 
iLebrun et al 12003) observed all sources except RX J0636. Hard 
(20-100 keV) X-ray data were extracted from all pointings 
within 12° from the source positions. 

To study the persistent X-ray emission, the time-averaged 
ISGRI spectra were obtained from mosaic images in five energy 
bands, logarithmically spaced between 20 and 100 keV. Data 
were reduced with the standard OSA software version 7.0 and 
then a nalyzed with the algorithms described bv lGoldwurm et aP 
(l200l . 



2.3. The Swift observations 

The Swift Burst Alert Telescope, BAT (lBarthelmvll2000l) . is a 
wide-field (-1 steradian) coded aperture mask instrument sensi- 
tive in the 14-195 keV range. Thanks to the large field of view, 
BAT has built up a sensitive all-sky map of the hard X-ray sky. 
The BAT team provides the average 8-channel spectrum of each 
source detected over the first 58 month of the mission. 

RX J0636 only recently appeared in the Swift/BAT 
58 month survey catalogue. To obtain an high- 
energy coverage of this source we downloaded 
the BAT spectrum along with the response matrix 
(http://swift.gsfc.nasa.gov/docs/swift/results/bs58mon/). 



3. Results 

3.1. Timing analysis 

All sources show periodic-like variability in their EPIC light 
curves and six of them also in their OM light curves. We then 
analyzed these light curves to identify periodic signals and to 
characterize them at diff'erent energies. 
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Table 1. Summary of main observations parameters for all instruments. 



Source 


Telescope 


OBSID 


Instrument* 


Date 
yyyy-mmm-dd 


UTstart 
hh:mm 


(ks)** 


Net count rate (mag)*** 
cts/s 


IGR J08390-4833 


XMM-Newton 
INTEGRAL 


0651540101 


EPIC-pn 
OM-B 
OM-UVM2 
IBIS/ISGRI 


2010 Dec 20 


12:34 
12:20 
16:57 


33.9 (27.3) 
15.0 
15.2 
2000 


0.776 ± 0.006 
3.6 ±0.9 (17.8) 
0.25 ± 0.05 (17.3) 
0.26 ± 0.04 


IGRJ18308-1232 


XMM-Newton 
INTEGRAL 


0601270501 


EPIC-pn 
IBIS/ISGRI 


2010 Mar 11 


15:26 


27.6 (15.7) 
2200 


0.793 ± 0.007 
0.36 ± 0.05 


IGR 116500-3307 


XMM-Newton 
INTEGRAL 


0601270401 


EPIC-pn 
OM-B« 
OM-UVWl" 
IBIS/ISGRI 


2010 Feb 2 


11:04 
10:50 
13:00 


30.0 (18.0) 
6.2 
8.0 
2200 


1.33 ±0.01 
10 ± 2 (16.7) 
2.2 ± 0.5 (16.3) 
0.31 ± 0.02 


IGR 118173-2509 


XMM-Newton 
INTEGRAL 


0601270301 


EPIC-pn 
OM-B 
IBIS/ISGRI 


2009 Sep 07 


02:03 
01:49 


35.0 (30.0) 
15.7 
4500 


0.941 ± 0.006 
3.0 ±1.0 (18.1) 
0.88 ± 0.04 


IGR J17195-4100 


XMM-Newton 
INTEGRAL 


0601270201 


EPIC-pn 
OM-B 
OM-UVM2'' 

RGS 1&2 
IBIS/ISGRI 


2009 Sep 03 


06:58 
06:44 
06:44 
06:35 


24.0° (27.0) 
11.4 
11.4 
33.9 
2600 


6.70 ± 0.02 
36 ± 5 (15.4) 
2.3 ± 0.6 (14.9) 

0.93 ± 0.05 


V2069 Cyg 


XMM-Newton 
INTEGRAL 


0601270101 


EPIC-pn 
OM-B 
OM-UVM2'' 
IBIS/ISGRI 


2009 Apr 30 


11:09 
10:55 
14:18 


26.4 (12.4) 
10.6 
10.6 
1900 


1.05±0.01 
10±1.5 (16.8) 
0.17 ±0.04 (17.6) 
0.14 ±0.05 


RX J0636+3535 


XMM-Newton 
Swift 


0551430601 


EPIC-pn 
OM-E" 
0M-UVM2« 
BAT 


2009 Mar 18 


17:01 
16:47 
22:41 


29.2 (12.1) 
6.7 
6.7 
-7500 


1.37 ± 0.01 
13 ± 2 (16.4) 
1.5 ±0.2 (15.4) 
0.00028 ± 0.00003 


IGR 115094-6649 


XMM-Newton 
INTEGRAL 


0551430301 


EPIC-pn 
OM-U 
OM-UVM2° 

RGS 1&2 
IBIS/ISGRI 


2009 Feb 02 


13:39 
13:25 
17:49 
13:16 


30.0 
14.2 
12.4 
31.9 
1400 


2.56±0.01 
36 ± 4 (14.3) 
2.8 ± 0.8 (14.7) 

0.61 ± 0.06 


XSS J0056+4548 


XMM-Newton 
INTEGRAL 


0501230301 


EPIC-pn 
OM-B 
OM-UVM2 

RGS 1&2 
IBIS/ISGRI 


2007 Dec 31 


04:56 
04:42 
07:01 
04:33 


15.0 (7.4) 
6.7 
6.7 
16.9 
240 


3.06 ± 0.02 
69 ± 5 (14.6) 
10 ±1(13.2) 

0.35 ± 0.06 



* Optical Monitor band reported when data usable 

** Net exposure times. In parenthesis is reported the solar flare removed exposure 
*** Instrumental magnitude reported in parenthesis. 
" No variability detected. 

* Time window filter applied: only first 2/3 of observation are used. 
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Table 2. Timing properties of the source sample. From left to right: (X-ray main (spin) period); P^-^^ (X-ray sidebands); Pq*"*^ (orbital period derived from X-ray sidebands); 
P^^^ (orbital period derived from phase-fitting); Af/A^.^^ (spin to sideband amplitude ratio); ^"f-'^'^ optical spin period from OM; P"^'^'^ (optical beat from OM); p^^'-^^'^'^ 
(optical spectroscopic period with its reference in parenthesis); (adopted orbital period in this work. 



Source 


pX 

w 










T}Opt,0M 


popt,OM 


nopt,spec. 


pA 




S 


S 


h 


h 




S 


S 


h 


h 


IGRJ0839 


1480.8±0.5 










1480±13 


1560±7 


8±1.0'' 


8±1 


IGRJ1830 


1820±2 


2400±19(w-2n) 


4.2+0.2 


- 


2.6 








4.2±0.2 


IGRJieSOO* 


571.9±0.5 














3.617+0.003 [1] 


3.617+0.003 


IGRJ1817 


831.7±0.7 (2w) 


810±1.0 (2(w + n)) 


8.5±0.2 


6.6±0.3 (3.41)^^ 


2.5 


1690±18 






6.3-8.7 


IGRJ1719'' 


1062±2 


1163±4(w-n) 


3.2+0.8 




1.4 




1102±4 


4.005±0.006 [1] 


4.005±0.006 


V2069 Cyg 


743.1+0.6 










751 ±5 




7.48039±0.00005 [2] 


7.48039+0.00005 


RXJ0636* 


920±1 


952±4 (£,; - Q.) 
1016+5 (6j-2n) 


7.6±0.9 
5.4±0.5 


10±i (3.00)'^'/ 


1.4 
1.9 






~ 3.4 [3] 


>5 


IGRJ1509« 


808.7±0.1 


842.4±0.6 (oj - n) 
412.7±0.5 (2a) - O) 


5.66+0.08 
5.6±0.3 




1.4 
2.4 


814±1 




5.89±0.01 [1] 


5.89±0.01 


XSSJ0056'' 


470.1 ±0.2 


454.2±0.7 (w + SI) 


3.7±0.2 
2.52+0.08' 


4.1±0.4(4.21)^ 


2.6 


474±2 




2.624±0.001 [4] 


2.624±0.001 



" Orbital period from OM photometric periods. 

* Ground based photometric optical pulse at 597.92+0.04 s [1]. 
^ Detection significance is reported in parenthesis in unity of cr. 

'' Ground based photometric optical pulse at 1 139.55+0.04 s and weaker signal at ~ 1055 s [1]. 

' Ground based photometric optical pulse at 1008.341 ±0.002 s and weaker signal at 930.583±0.004 s [3]. 

^ Since the total exposure is ~8 h, we here give the uncertainty at 3cr c.l. 

* Beat period from 3(6; - Q). Amplitude ratio refers to 3a) and 3(a) - f2). Ground based photometric optical pulse at 809.42±0.02 s [1] 

* OM B band pulses derived from third and fifth segment. Ground based photometric optical pulse at 470.40±0.03 s [5]. 

' Orbital period from X-ray light curve. 



References: [1] Pretorius 2009; [2] Thorstensen et al. 2001; [3] Gansicke et al. 2005; [4] Bonnet-Bidaud et al. 2009; [5] de Martino et al. (privcomm.) 
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XSS J0056+4548 



0.5 



0.4 - 



P=2.6(4) h 




0.2 



0.1 



0.2 



0.25 0.3 0.35 

Time since BJD=2454465.5 [daysj 



Fig. 4. PF vs time for XSS J0056. The sinusoid period with a 3cr 
uncertainty is also reported. 



For each source, the Hght curve in the 0.3-15 keV range was 
folded at the main (P„) period reported in table |2] and then fitted 
with a series of sinusoids (one for each harmonic). The statistical 
significance of the inclusion of higher harmonics with respect 
to the fundamental one was evaluated by an F-test (see table 
[3]l. The first harmonic was always statistically significant (here- 
after we consider significant a probability greater then 3cr c.l.), 
with the exception of IGRJ1830 and IGR J1719 in which only 
the fundamental is significant. The second harmonic (3co) was 
found to be significant in IGR J0839, V2069 Cyg and IGR J1509. 
V2069 Cyg is also the only source showing 4a>. 

The power spectra of most sources either show an excess of 
power around the fundamental frequency or additional weaker 
peaks close-by. In some cases also the harmonics show similar 
excess or peaks. Therefore, to estimate their frequency we de- 
trended the light curves from the main pulse by subtracting the 
composite sinusoidal function corresponding to the fundamen- 
tal frequency plus their harmonics when detected. The residual 
power spectra of the detrended light curves makes clearly visi- 
ble additional peaks. Consequently, we fitted the detrended light 
curves with one new sinusoid (or more sinusoids when harmon- 
ics are present). We reported the estimated period, when they re- 
sulted to be statistically significant, in table |2] These additional 
periodicities are identified as the orbital sidebands of the main 
pulse (lo). 

IGR J 17 19, RXJ0636 show the negative sideband w - O 
(beat). In IGRJ1509 we detected the harmonic of the beat 
3(L) - Q, but not the fundamental, and the 2a) - Q sideband. In 
RXJ0636 we also detected oj - 20. For this source the pres- 
ence of both negative sidebands is significant at 4cr c.l.. For 
IGRJ1830 the w - 2Q is detected at 4.1cr c.l., but the w - O 
is only at 1.5cr c.l.. The positive ai + Q. sideband is detected in 
XSS J0056 (7.4o- c.l.) and 2(w -I- O) in IGR J 1 8 1 7 ( 1 6 o-c.l.) since 
the main X-ray period is the first harmon ic (see sect. 13.1.3) . 

We here note that ISazonov et al. I (2008) using a short 
Chandra observation of IGRJ0839 detected a periodicity of 
1450 + 40 s consiste nt with our more constrained value. 
iButters et al] (|2008|) failed to detect a periodicity in a RXTE 
observation of IGRJ1719. On the other hand, they detected 



two periodic signals in XSS J0056 again using RXTE. One at 
465.68+0.07 s and a weaker one at 489.0+0.7 s. The former is 
inconsistent with our determination (being out by more than 5 
s), while the la tter is not detected i n our data. 
For IGR J 1 8 H lNichelh etaP (12009) using a S wift/XRT pointing 
report a period of 830.7 ±0.07 s consistent with our results and a 
period of ~ 1660 s from a previous short Chandra exposure (see 
sec t.l3.1.3]for its i nterpr etation). 

IButters et all (l2009l) by using RXTE data of IGR J 1509 
found a signal at 809.7+0.6 s, consistent with our results. 

The phase-fitting procedure also allowed us to get an esti- 
mate of the orbital period. The light curve of each source was 
subdivided in n intervals (10, 15 or 20), where n depends on 
the source brightness and on the intensity of the signal. Then 
we folded each light curve at the period reported in table 
|2] and then we fitted it with a Fourier sine series truncated at 
the highest significant harmonic. We consequently studied the 
variations with time of the phase (and amplitude) of the main 
pulsation, corresponding to the fundamental frequency. In most 
cases the phase variations, typically within A<p+ 0.2, appear to 
be randomly distributed around the average value. Only in three 
sources, IGR J1817, RX J0636 and XSS J0056 the lag variations 
show a clear trend, allowing a sinusoidal fit. These are shown in 
figure[Tl|2] and[3] The statistical significance of the inclusion of 
a sinusoid respect to a constant function is evaluated with an F- 
test, and was found to be 3.41, 3.00 and 4. 21cr respectively. The 
resulting periods are reported in column 5 of table |2] 

Most sources show also variations with time of the pulse am- 
plitude (or PF, see sect. 13.1.21 for its definition), but these ap- 
pear to be randomly distributed around an average value, with- 
out showing a clear trend. An exception is XSS J0056, where 
also the amplitude changes following a sinusoidal modulation. 
A sinusoidal fit is statistically significant at 3.52cr (see figure|4|i, 
giving a period of 2.6 + 0.4 h (3crc.l.). Hence, while in RX J0636 
and IGR J1817 the phase shifts could be due to time arrival de- 
lays, in XSS J0056 the variations seen in both amplitudes and 
phases are likely intrisic. This is also suggested by the presence 
of a peak around 9.3 x 10"^ Hz in the power spectrum and a si- 
nusoidal fit gives a period of 2 .52+0.08 h, in agreement wi th the 
spectroscopic period found bv lBonnet-Bidaud et al. I (l2009h . 

Interpreting the main pulsation as the WD spin period, the 
sidebands can be used to estimate the binary orbital periods, re- 
ported in column 4 of table |2] For IGR J 17 19 and IGR J 1509 
the periods are in agreement with those derived from spec- 
troscopy, reported for comparison in column 9. The values found 
for IGR J1817 and XSS J0056 are consistent within errors with 
those derived from the phase-fitting method, but for XSS J0056 
they are not compatible with the results from the PF vs time anal- 
ysis and from the sinusoidal fit. For RX J0636 we found con- 
sistency between the period found from phase-fitting and that 
derived using the beat but not a> - 2f2. However, all derived 
values are largely off from the tentative period of 3.4 h found 
from optical spectroscopy ( Gansicke et al.ll2005 ). Furthermore, 
V2069 Cyg also shows in its power spectrum substantial power 
at low frequencies, but we are unable to constrain a period from 
the light curve. However, a sinusoidal fit with period fixed at 7.5 
h, derived from spectroscopy by Thorstensen & Tavlor (200l|, 
provides a possible match with data and gives a 35% amplitude 
modulation. This variability is more pronounced at soft energies 
(<3keV), where we also detect a dip, lasting ~1500 s, around 
the minimum of this long-term X-ray modulation (figure |5). 

For the systems under study we hereafter adopt, when avail- 
able, the more accurate orbital spectroscopic periods, except for 
RXJ0636 for which the orbital period is likely larger, Pn ^5 
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Fig. 1. Left Panel: Power spectrum. Fundamental and harmonics are reported in red, while dotted blue lines represent sidebands. 
Right panel: Evolution of the phase of the main signal, w, vs time. 
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Table 3. Pulsed fraction vs energy. Highest detected harmonic and its statistical significance is also reported. 



Source 






Pulsed Fraction 




na> 


sign. 




0.3-1 keV 


1-3 keV 


3-5 keV 5-15 keV 


0.3-15 keV 








% 


% 


% % 


% 




cr 


IGR J0839 


99 ±2 


74 ± 1 


35 ±2 18 ±2 


57 ± 1 


3 


3.97 


IGR J1830 


30 ±4 


33 ± 2 


26 ±2 21 ±2 


28 ± 1 


1 




IGR J1650 


18 ±2 


12 ± 1 


5±2 6±2 


9± 1 


2 


4.29 


IGR 11817" 


58 ±3 


46 ±2 


34 ± 2 27 ± 1 


30 ± 1 


2 


5.94 


IGR 11719* 


14 ± 1 


13 ± 1 


6± 1 < 3 


10 ± 1 


1 




V2069 Cyg 


14 ±3 


25+2 


16 ±1 9±1 


18 ± 1 


4 


3.5 


RX J0636 


18±2 


21 ±2 


< 17 


15 ± 1 


2 


3.02 


IGR 11509 


57 ± 1 


51 ± 1 


27 ± 1 9 ± 1 


42 ± 1 


3 


3.52 


XSS J0056 


42 ±2 


31 ± 1 


16 ±2 9±1 


24 ± 1 


2 


4.54 



" Pulsed fraction refers to 2cj. 

* Derived from the first 2/3 of observation. 
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Fig. 5. 0.3-1 keV hghtcurve of V2069 Cyg binned at 198 s. 
Reference time is the PN start of exposure, see table [T] A dip 
is present at about 1 .5 x 10"* s. 



h. For those sources without a spectroscopic determination we 
adopt the orbital periods as derived from the X-ray sidebands 
(for IGR J0839 from optical sideband). Exception is IGR JIB 17, 
for which we adopt the range 6.3-8.7 h, resulting from sidebands 
and phase-fitting method. The adopted values are reported in the 
last column of tabled 

3.1.2. Pulsed fraction 

We also studied the energy dependence of pulses. We selected 
five energy bands, producing four background subtracted light 
curves binned at 66 s (0.3-1, 1-3, 3-5, 5-15 keV), plus the 0.3- 
15 keV light curve with a bin time of 11 s, and folded them 
at the value of P^ reported in table |2] For all sources the X- 
ray pulse profiles have structured pulse shapes. Moreover, for all 
sources the pulse profiles resulted to be roughly phase aligned at 
all energies, within statistical uncertainties (see figure |6] and [T). 



The folded light curves were then fitted with a Fourier sine 
series truncated at the highest significant harmonic. This allowed 
to study, for each energy interval, the pulsed fraction (PF), here 
defined as: PF = (Amax - Ainin)/(Ainax + Kmi), where A^ax and 
Amin are respectively the maximum and minimum value of the 
sinusoid used to account for the fundamental harmonic. For all 
sources the PF was found to clearly decrease with energy (see 
table|3]l. In order to further inspect for possible spectral changes 
with the spin phase, we also construct the Hardness Ratios (HRs) 
using the three energy intervals 0.3-1 keV 1-3 and 3-5 keV. 
The HR is here defined as: HR = iij,gy;,/n]^^y^, where nj^^^ is the 
number of photons detected by the PN camera in the inspected 
energy range (called a and b), and / is the phase interval. From 
HR = Uj^/Uj 3 we found for all sources a hardening at the 
pulse minimum, except for IGR J 1830 and V2069 Cyg, where 
HR was found to be constant. For V2069 Cyg, due to its peculiar 
spectrum in the soft energy range (E< 3 keV; see sect. 13.2. fl i. we 
also inspected the HR = n'j j/Uq ^ j, that indicate a softening at 
spin minimum, suggesting the contribution of a soft component. 



3.1 .3. The Optical and UV variability 

Optical and UV light curves were inspected for the presence of 
periodic signals by producing power spectra in each band, except 
for RX J0636, IGR J 1 830 and IGR J 1 650 that were too faint. The 
UV light curves are generally noisier and for most of them we 
are unable to find significant periodicities. 

The light curves were also fitted with a sinusoidal function. 
The periods were always found to be consistent with one of those 
detected in the X-rays (see table|2]l. Worth noticing are the fol- 
lowing cases: 

IGR J0839: The UV Hght is modulated (PF=35%) at a period 
(1480 + 13 s) consistent with that found in the X-rays, but the 
B band light curve shows a period of 1560 + 7 s (PF=22%). 
The latter is likely the orbital sideband (beat) w-Q, suggesting an 
orbital period of Pn=8 + 1 h. Hence this source is spin-dominated 
in the X-rays and UV, but beat-dominated in the optical. 

IGR J1817: From the light curve in the B filter we derived a 
period of 1690+18 s, twice the detected X-ray period. Moreover, 
no signal was found at the X-ray period. We conclude that the 
X-ray period corresponds to the first harmonic, while the optical 
period is the spin. 
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Fig. 6. X-ray pulse profiles at different energy intervals for different sources. Energy increases from top to bottom. Optical and UV 
band pulse profiles are also reported when available. 
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Fig. 8. OM light curve of XSS J0056 in the B band for the five 
time intervals of exposure. Time increases from bottom to top. 



IGRJ1719: The B band light curve, simultaneous with the 
X-rays, shows a period of 1102 + 4 s. This value falls in be- 
tween the two periods detected in the X-rays (1062 + 2 s and 
1 163 + 4 s), suggesting that both periodicities are also present at 
optical wavelengt hs. This is also c onfirmed by ground-based op- 
tical photometry (|Pretoriusll2009l) that showed a weaker ~ 1055 
s and a stronger ~ 1 139 s peak in the power spectra. Therefore, 
also this system is spin-dominated in the X-ray, while beat- 
dominated at optical wavelengths (see also sect. 14.21 ). 

V2069 Cyg: we derived a pulse period of 751 + 5 s, consis- 
tent, within errors, with that found in the X-rays. Interestingly, 
folding the light curve at the X-ray period, the optical pulse is 
anti-phased with the X-ray one. Moreover, it shows a broad max- 
imum centered on the X-ray minimum, with a dip where an X- 
ray secondary maximum occurs (see Fig|6). 

XSS J0056: it has peculiar UV and optical variabilities. The 
power spectra in these bands did not show significant peaks. 
However, inspection of the individual five segments (~1.34ks 
each) in the two bands revealed the presence of non-stationary 
pulses. They are found to be phase aligned with the X-ray pulse 
during the first 2/5 of the OM B band exposure, for which a sinu- 
soidal fit gives a period of 474+4s (consistent with what found 
in X-rays), but antiphased in the third and fifth segment (see fig- 
ure Hand [8]i. The UVM2 filter data revealed a weak variability 
(~ 9%) during the third and possibly the fourth segment, but we 
are unable to determine the period. 

3.2. Spectral analysis 

3.2.1 . 0.3-1 00 keV spectrum 

For all sources, a single-temperature model is inadequate to de- 
scribe the broad band (0.3-100 keV) spectrum. We found the 
best match with composite models that consist of a combination 
of a number of optically thin plasma (mekal) components plus a 
Gaussian at 6.4 keV, including a total (Wabs) and a partial cov- 
ering absorber (Pcfabs). The statistical significance of the inclu- 
sion of additional mekal components was always verified with a 
F-test. All best-fit spectral parameters are reported in table|4] 
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At least two characteristic temperatures were always 
present, three in the case of IGRJ0839, IGRJ1830, IGRJ1650, 
IGRJ1817 and IGRJ1719. We broadly divided these optically 
thin components into a cold (kTc ~ 0.1 - 0.2 keV, exception is 
IGR J0839 with kTc ~ 0.7 keV), a medium (kT„ -4-10 keV) 
and an hot (kTh ~30-60keV) one. The presence of low tem- 
perature plasma is supported by the RGS spectral analysis of 
IGR J1719, IGRJ1509, and XSS J0056 (see sect. [TO] ). 

However, in V2069 Cyg the cold mekal temperature was un- 
constrained, while for RX J0636, the abundance was unreliably 
low (Az ~0) A more physically-justified model is obtained by 
substituting the cold mekal with a blackbody (BB) component 
at a temperature of kTsB ~ 70-80 eV (see also figure|9]l. 

All metal abundances were found to be consistent with the 
solar value except IGRJ1817, for which a highly sub-solar 
(0.20 + 0.07) value was found. 

For all sources we found significant presence of an emission 
feature at about 6.4 keV (see figure |9), with no significant shift 
in energy. Hence, we fixed the centroid energy of the Gaussian 
component at 6.4 keV. The equivalent widths (EW) of this emis- 
sion feature were found ~1 00-220 eV, consistent with fluores- 
cent Fe K„ line. This suggests the presence of a continuum re- 
flected component, whose inclusion (reflect), however, pro- 
vides meaningful fit only for IGR J 1650, lowering the tempera- 
ture of the hot MEKAL (see last line of tablelDl. 

All sources are affected by strong complex absorption. The 
hydrogen column density (Nh) of the total absorber was found 
to be lower or consistent with the interstellar value in the di- 
rection of the sources. The partial absorber, with typical cov- 
ering fraction ~ 60%, was found at much higher column den- 
sities (up to lO^-'cm"^) and hence is instead intrinsic to the 
sources. For IGR J1817 and XSS J0056, which were found with 
the highest value of intrinsic absorption (50 + 10 x 10^^ and 
300 ± 200 X 10^^ cm'^) we used the Pwab XSPEC model 
dPone & Magdziarz 1998), that assumes a power-law distribu- 
tion of neutral absorbers and is likely more appropriate in case 
of high density intervening absorbing columns. This provided a 
better fit in terms ofx^. However, the column density resulted to 
be poorly constrained in XSS J0056. 

IGR J0839 is the only source which, in addition to the two in- 
tervening absorbing columns, also displays an evident absorp- 
tion feature found at 0.76 keV. We modeled it using an edge 
component, which statistical significance, evaluated with an 
F-Test, resulted to be 5.6cr. This feature is consistent with being 
an OVII K-shell absorption edge with an optical depth t=0.6. It 
indicates an ionized absorber along the line of sight. 
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Fig. 9. Broad-band unfolded spectra, with post fit residuals in the 
lower panels, for three representative cases. Spectra have been 
rebinned for plotting purpose. Composite model (solid black 
line), individual components (dotted blue lines). Up: IGR J1719. 
Wabs*Pcfabs*(3mekal-i-gauss) model used. Central: V2069 Cyg. 
Wabs*Pcfabs*(bbodyrad-i-mekalH-gauss) model used. Bottom: 
IGR J0839. Wabs*Pcfabs*(3mekal+gauss) model used. An ab- 
sorption feature at ~0.76 keV is present in the residuals (lower 
panel). 
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Table 4. Spectral properties for the best fitting model. 
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2.04(2) 


2.3(2) 


1.23/900 






























XSS J0056 


2mek* 


0.04(1) 


300(200)^ 


-0.60(1)'' 


0.16(1) 


10.2(6) 




9.0(4) 


50(10) 




0.42(5) 


0.13(3) 


3.4(3) 


1.8(4) 


^Vdof 1.03/578 






























IGRJ1650 


R*3mek« 


0.28(2) 


9(1) 


60(1) 


0.18(2) 


7(1) 


41(4) 


0.4(2) 


1.1(2) 


6(3) 


1.6(fix) 


0.22(2) 


1.5(1) 


2.8(2) 


;i^2/dof 1.08/630 































" E,i,, = 0.75(l),Ted, = 1.3(2). 

* Pwab model used. Nh^,„ is fixed at the minimum boundary level of 10'^ cm ^. 
" /3. 

" Blackbodyrad kT. 

^ Blackbodyrad normalization xlO**. It is defined as n = R^/Djp, where Dio is the source distance in unit of 10 kpc. 
« Reflection model used. R=0.98±0.02. Cosine of incUnation is fixed at 0.45. 
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Table 5. column 2: WD m asses as derived fr o m ma ximum post- 
shock temperature, using ISuleimanov et al.l (l2005b model (see 
text), column 3: Lower limits to the mass accretion rate adopting 
the source distances reported in column 4 and the orbital periods 
in the last column oftable\2\ 



Source 


MwD 


M 




D* 




Mo 


Mo yr 


1 


pc 


IGR J0839 


0.95+0.08 


3.2x10- 


10 


1200 


IGR J1830 


0.85+0.06 


2.4x10" 


-11 


100 


IGR J 1650 


0.92+0.06 


4.1x10" 


-11 


270 


IGR J1817 


0.96+0.05 


2.4x10" 


10 


610 


IGRJ1719 


0.86+0.06 


1.9x10" 


11 


130 


V2069 Cyg 


0.82+0.08 


5.3x10" 


-10 


930 


RX 10636 


0.74+0.06 


5.4x10" 


-10 


1100 


IRGJ1509 


0.89+0.08 


2.5x10" 


-10 


600 


XSS J0056 


0.79+0.07 


2.5x10" 


-10 


320 



* Lower limits to the distances are determined as follows: we used the 
lower uncertainty values of the adopted orbital periods and the 2MASS 
K band (reddening correction not applied) magnitudes of the optical 
counterparts. For the secondary stars we used the K band absolute mag- 
nitudes predicted for CV donors as a function of P^rb (Knigge 2006), 
assuming that these donors totally contribute to the K band flux. 



Since the hot temperature should be regarded as a lower limit 
to the maximum temperatu re of the PSR, we also used the model 
of lSuleimanov et al](l2005h . that takes into account the growth of 
pressure toward the WD surface and hence the change of gravity, 
to obtain a more reliable value for the maximum temperature and 
consequently an estimate of the WD mass. This model is com- 
puted for the continuum only. We included in the fit the same 
absorber components as in the multi-temperature fit, as well as 
a broad Gaussian that takes into account the iron line complex. 
We obtain in all cases an acceptable between 1.1-1.3. The 
resulting masses are found in the range 0.74-0.96 Mq (see table 
|5]l. We also note that the temperature of the hot component (and 
hence the masses) determined with the multi-temperature model 
are within 3cr consistent with those found with the PSR model, 
except for XSS J0056 whose temperature would give a unreal- 
istic low WD mass. This could suggest that the hot component 
is close to that at the shock. All other best fit parameters are 
compatible, within uncertainty, with that reported in table |4] For 
IGR J 1650, we also tried to include a reflection component with- 
out any fit improvement. This is likely due to the lack of cover- 
age in the 10-20keV range and the low statistics above 20keV, 
that prevent us from drawing conclusions on the reflected con- 
tinuum. 

3.2.2. Pulse Phase Spectroscopy 

To understand the role of spectral parameters in producing the 
observed X-ray pulses, we carried out a phase resolved spectral 
analysis. For all sources we extracted the spectra at the maxi- 
mum and minimum of the pulsation (see table |6] for details on 
the interval selection). We kept fixed the column density of the 
total absorber, the metal abundance, and the parameters of the 
hot MEKAL component to the values found from the phase av- 
eraged spectral fits. We also fixed the intermediate mekal tem- 
perature, since it was always consistent with that derived from 
average fit (except for IGRJ1650). Furthermore, inspection of 
the fluorescent Fe line for phase-dependent velocity shifts gave 
no significant change and hence we fixed the Gaussian centroid 



at 6.4 keV. All other parameters were left free to vary (see table 

In all cases the spectrum at the minimum is more absorbed 
than that at the maximum (see table |6]l, due to higher values of 
the hydrogen column density and of the covering fraction (how- 
ever for XSS J0056 we fixed the density of the absorber to the 
average spectrum values). No variation in the normalization of 
the cold MEKAL component is found except for IGR J 1650 where 
it decreases at minimum. We then conclude that the pulses are 
mainly due to phase-dependent absorption. 

Interesting is the case of V2069 Cyg that instead shows an 
increase in the normalization of the BB component at pulse min- 
imum, in accordance with the behaviour of the energy dependent 
pulses (figure |6ll. This suggests that the secondary maximum at 
(p ~ 0.8, best seen below 1 keV, is produced by an increase of the 
normalization of this optically thick component. 

We found evidence for phase variation in the EWs of the 
Fe fluorescent fine only in IGRJ0839 and IGRJ1509. In these 
sources the EWs change between 0.08+0.03 and 0. 16+0.03 keV 
and 0.09 + 0.03 and 0.23 ± 0.04 keV at maximum and minimum, 
respectively. 

3.2.3. RGS spectrum 

The RGS spectra of IGRJ1719, IGRJ1509 and XSS J0056 (see 
figure |9j were fitted in two ways. In one, we used continuum 
only models (absorbed BB plus bremsstrahlung) and fitted the 
strongest lines using Gaussians, to derive their centroid, width, 
and intensity in a manner that does not strongly depend on 
the choice of models. In IGR J1719, the fines of OVII, OVIII, 
NelX and NeX were all strong enough for such a treatment; in 
IGR J1509, only the two oxygen fines; and for XSS J0056 only 
the OVII line (see table |7]). In addition to the line fluxes, which 
are consistent with the prediction of the EPIC model, we found, 
for IGRJ1719 and IGR J 1509, that the fines are mildly broad- 
ened with Gaussian cr of order 1000 km s i (or FWHM -2500 
km s ') and that the centroids of the He-like triplets (OVII and 
NelX for the former and OVII only for the latter) are between the 
expected energies of the resonant and intercombination compo- 
nents, when fitted with a single Gaussian. 

We also explored global fits, based on the best-fit model for 
the EPIC spectra. For XSSJ0056, there is a marginal indica- 
tion that an addition of a OVII edge would improve the fit, even 
though this was not statistically significant in the EPIC fit. This 
should be explored further using future, higher signal-to-noise 
spectra. For IGR J 17 19, although the basic features of the RGS 
spectra were reproduced using the EPIC model (including the 
normalization), the reduced is 1.32 for 363 degrees of free- 
dom. We obtained an improved fit (x^ - 1.03, 351 dof) by in- 
cluding a Gaussian deconvolution using the XSPEC g smooth 
model and by letting the mekal density parameter be fitted. The 
fit converged to a high density regime (i.e., with no forbidden 
components for the He-like triplets) with cr ~ 1250 km s 
consistent with our findings from the phenomenological model. 
Similarly, a high density, mildly broadened model was more suc- 
cessful in fitting the RGS spectra of IGR J1509. 



4. Discussion 

The CVs analyzed here are all found to be strong X-ray pulsators 
with periods ranging from -470 s to -1800 s. In IGRJ0839, 
IGRJ1817, IGR J1719, V2069Cyg, IGR J1509 and XSS J0056, 
periodicities are also identified in simultaneous OM optical (B 
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Table 6. Spectral parameters of max and min spectra. Other parameters are fixed to their best-fit values of the average spectrum. 



source 


10^^ cm"^ 


cvf 

% 


kTe 
keV 


He 

10"^ 


lO""* 


10"' 


EW 
keV 


Fo.3-15 

10"" 




IGR 10839° 
Max (0 = 0.0 - 0.4) 
Min (0 = 0.5 - 0.9) 


6(1) 
7(1) 


40(3) 
81(1) 


0.69(fix) 
0.69(fix) 


0.34(7) 
< 0.4 


1.2(3) 
0.7(3) 


3.9(3) 
3.1(2) 


0.08(3) 
0.16(3) 


0.92(3) 
0.60(3) 


1.20/485 


IGR 11830 
Max (0 = 0.9 - 1.1) 
Min (0 = 0.35 - 0.65) 


5(2) 
14(2) 


49(3) 
65(2) 


0.21(1) 
0.16(1) 


5(1) 
30(10) 


1.9(3) 
1.3(2) 


4.7(3) 
3.9(5) 


0.22(4) 
0.23(8) 


1.11(6) 
0.76(4) 


1 m / 1 T 1 
l.UJ/l 31 


IGR 11650 
Max (0 = 0.8 - 1.1) 
Min (0 = 0.2 - 0.5) 


13(2) 
12(1) 


56(2) 
66(1) 


0.09(1) ; 0.7(1)° 
0.15(4); 8(2)'' 


2.8(4) 
0.6(2) 


0.10(5) 
1.8(5) 


8.6(5) 
6.6(7) 


0.25(5) 
0.22(3) 


1.58(4) 
1.44(6) 


0.93/353 


IGR 11817 
Max (0 = 0.55 - 0.75) 
Min (0 = 0.0 - 0.2) 


50(15)^ 
80(20)'- 


-0.25(8) 
-O.Ul)"* 


0.26(fix) 
0.26(fix) 


9.6(3) 
16(7) 


15.2(6) 
13(9) 


10(4) 
9(3) 


0.15(4) 
0.10(3) 


1.7(1) 
1.4(1) 


1.06/356 


IGR 11719 
Max (0 = 0.2 - 0.5) 
Min (0 = 0.7-1.0) 


6.7(7) 
7.9(6) 


39(2) 
47(1) 


0.18(1) 
0.15(1) 


1.2(1) 
1.7(4) 


6(1) 
5(1) 


19(1) 
20(1) 


0.14(2) 
0.13(2) 


4.75(3) 
4.64(3) 


1.03/1768 


V2069 Cyg 
Max (0 = 0.1 - 0.55) 
Min (0 = 0.6 - 0.95) 


5.4(6) 
11(1) 


61(7) 
70(1) 


0.073(1) 
0.065(1)^ 


5.8(6)' 
13(1)^ 


- 


7.2(1) 
6.6(1) 


0.24(8) 
0.27(4) 


1.33(3) 
1.07(4) 


1.05/275 


RX J0636 
Max ("rti - 3 - 7'! 
Min (0 = 0.8-1.1) 


9('21 
8(2) 


64(3) 


08ny 

0.08(1)" 


6r 1 \f 

0.7(3)' 




4.7(3) 


osn'i 

0.07(6) 


1.04(3) 
0.87(5) 


1.08/441 


IGR 11509 
Max (0 = 0.75- 1.1) 
Min (0 = 0.2 - 0.6) 


5(1) 
9(1) 


27(2) 
67(1) 


0.11(1) 
0.12(1) 


3.2(6) 
2.8(6) 




11.6(3) 
10.7(3) 


0.09(3) 
0.23(4) 


2.50(4) 
1.82(2) 


1.11/824 


XSS J0056 
Max (0 = 0.4 - 0.7) 
Min (0 = 0.0 - 0.2) 


300(fix)'- 
300(fix)'- 


-0.64(1)'' 
-0.52(1)'' 


0.16(2) 
0.12(4) 


8(2) 
20(10) 


50.0(2) 
48.6(2) 




0.12(3) 
0.15(3) 


3.2(1) 
2.6(1) 


0.96/360 



" Max: E^rfj, = 0.75 ± 0.01 r = 1.4 + 0.2; Min: Erfg^ = 0.76 ± 0.01 t = 0.9 ± 0.2. 
* Medium mekal temperature. 



' Blackbodyrad temperature. 

^ Blackbodyrad Normalization xlO''. Blackbodyrad normalization is defined as R^^^/Dj^, where Dio is the source distance in unit of 10 kpc. 



or U bands) or UV photometry. X-ray pulses are unambiguous 
signature of magnetically confined accretion and hence of the 
WD rotational period. Therefore, all sources can be classified as 
CVs of the IP type. 

The XMM -Newton data have allowed us: to identify for 
the first time the WD spin period in IGRJ1830, IGRJ1650, 
IGR Jl 8 17, IGR J 17 19, V2069Cyg and RX J0636, and to refine 
the spin period of IGR J0839 and IGR J 1509 (but we detemiined 
a slightly different period for XSS J0056). Furthermore, we 
could infer the orbital periods (except for IGR J1650). We pro- 
vide first estimates for IGRJ0839, IGRJ1830, IGRJ1817 and 
find consistency with the spectroscopic periods for IGR J 17 19, 
IGR J1509, V2069Cyg and XSS J0056 except RX J0636. 

We have also presented the first broad band X-ray spectral 
analysis of these sources as well as spin phase-resolved spectra, 
allowing us to identify commonalities to IPs, but also peculiari- 
ties in some sources. 

We discuss here the results in terms of accretion and emis- 
sion properties. 

4. 1 . X-ray pulses as tracers of accretion mode 

The presence in the power spectra of a strong and dominant sig- 
nal at the WD spin period, especially in the hard (E>4 keV) X- 



rays, is a signature that accretion occurs predominantly through a 
disc bec ause the circulating m aterial loses memory of the orbital 
motion ( Wvnn&King|ll992 '). For IGRJ0839, V2069Cyg and 
IGR J 1 6500, we have detected only the spin frequency in X-rays, 
and then they likely are pure disc-accretors. Noteworthy, the OM 
optical (but not the UV) data of IGR J0839 and the g round-based 
optical photo metry of IGR J 1 6500, IGR J 1 7 1 9 (Pre toriusll2009l) 
and RX J0636 toansic ke et aDl2005h reveal instead stronger sig- 
nals at the orbital sideband, w - Q, rather than at the spin period. 
Since the optical light originates from reproces sing of X-rays at 
fixed regions within the binary (lWarnerl[l995l) . this makes the 
X-rays the best tracers of the true WD rotational period. 

Six out of nine systems, while displaying a dominant spin 
periodicity in the X-ray band, also show signals at sideband fre- 
quencies (table|2]), with spin-to-sideband amplitude ratios ~ 1 .4- 
2.6. Hence, although the bulk of accretion onto the WD takes 
place through a disc, a fraction of it (~35% up to ~70%) passes 
over the disc and is directly channeled onto the WD magnetic 
poles. This hybrid mode of accretion, the disc -overflow, is also 
observed in other confir med IPs (Hellier 1 995tlNorton & Tavlod 
[1996'; Norton et al."l997) and in more recently identified mem- 
bers (Anzolin et al. 2008). Worth mentioning is the case of 
IGR J 17 19 whose X-ray light curve reveals spin and sideband 
pulses only during the first 2/3 of observation, while no strong 
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Table 7. Spectral parameters with RGS. 





1GR11719 




Line 


E 


tr 


Flux 




keV 


eV 


photons cm"^ s"' 


OVll 
OVllI 
NelX 
NeX 


5706+° '"" 

"■-"""-0.0008 

6516+"'""* 

"•"-"^"-0.0013 

0.9209+S 
1.0228±°:S 


3.8 + 0.7 

^•■^-1.4 
11 3+** 5 

< 32.2 


(3.3+^5) X 10-3 
(7.7±]|) X lO-'' 
(3.7±|^'*) X lO"'* 
(1.3±+'8") X 10 * 


o"s = 


1250±^^^ [km/s] 








ICR J 1509 




Line 


E 


cr 


Flux 




keV 


eV 


photons cm"^ s"' 


OVll 


5705+'''"'2i 

"■-" "-^-00015 


4.1±j;3 


(7.3±5 *) X 10"' 


OVllI 


0-6549+™ 


2.8±t? 


(3.6±i|) X 10"' 


lT„ = 

s 


1240±^;;« [knVs] 








XSS0056 




Line 


E 


(T 


Flux 




keV 


eV 


photons cm"^ s"' 


OVII 


5676+""*"* 

"■-'"'"-0.0444 


< 240 


(3.7±J"/) X 10-* 



coherent signals are found afterward. This system is the one with 
a low spin to beat amplitude ratio of ~1.4 and a spin PF~ 10%. 
It could be possible that the two variabilities interfere destruc- 
tively, as also suggested by the large jumps in the phase of the 
signal with time (see figure|2]l and the low-frequency variations 
in the power spectrum. Hence, the somewhat peculiar behaviour 
of this IP could be due to substantial disc-overflow. 

Two systems, IGR J1817 and XSS J0056, reveal a weak sig- 
nal at a putative positive sideband, to + Q., that would imply a 
retrograde WD rotation. We cannot exclude that the true spin 
period is the latter and hence the stronger signal is the beat. In 
this case the two systems would be beat-dominated at both X- 
ray and optical wavelengths and then disc-less accretor. A simi- 
lar problem raised for BG CM i (iNorton et al.ll992i iGarlick et all 
[19941 Id^ Martino et al To solve this issue, both fast time- 

resolved spectroscopy and polarimetry are needed. 

Furthermore, for IGR J 1830 and RXJ0636 the presence of 
the to - 2Q sideband would imply an orbital modulation of the 
amplitude of the beat frequency if this is an intrinsic emission. 
Hence, the beat and the spin pulses would be produced by dif- 
fere nt modes, stream-fed and disc-fed accretion respectively, see 
also [Norton et all (Il997h . 

Long-term variation consistent with the orbital period was 
found in XSS J0056, and likely in V2069Cyg. In the latter, we 
also detected a dip at minimum of the possible orbital modu- 
lation, more pronounced in the soft bands, indicating a spectral 
hardening. Orbital variations are often found in IPs and are likely 
due to absorption effects from fixed regions within the binary 
frame, such as the disc edge, in systems viewed at relatively large 
inclination angles, i^ 60" (Parker et al. 2005). 

Using the adopted orbital periods, reported in table |2] the 
spin-orbit period ratios, defined as P^/Pn, are lower than 0.1, 



(a) 




Channel Energy [keV] 



Fig. 10. From top to bottom, IGRJ1719, IGRJ1509, and 
XSS J0056. RGS 1 data and model (black), RGS2 data and model 
(blue). The models use the EPIC best-fit values. Those for 
IGR J 17 19 and IGR J 1509 also include fitting for g smooth cr 
and for mekal density parameters. Top panel inset shows the 
OVII line wavelength range (20.6-23.1 A) with the EPIC best- 
fit low density mekal (black), the high density, narrow mekal 
(blue) and the high density, gsmooth broadened meakl (red). For 
XSS J0056 we show the EPIC best-fit model with (black) and 
without (red) the edge. 



except for IGRJ1830 (0.12), and cluster around 0.03-0.05 (see 
figure[TTll- All but one (XSS J0056), have Pn above the orbital 
gap. A strong asynchronism (P^^/Pn < 0.1) and long orbital pe- 
riods are predicted by evolutionary models, where mCVs evolve 
towards shorter orbital periods and, depending on the WD mag- 
netic field strength, may or may not eventually synchronize 
(iNorton et al.ll2004l l2008h . In long period highly asynchronous 
systems, the accretion flow is expected to take part in a disc-like 
configuration, as indeed inferred from X-rays. Such configura- 
tion may change while the system evolves towards short orbital 
periods, depending on the WD magnetic field strength. To fur- 
ther progress in this issue, polarimetric measures of new systems 
are crucial to assess whether these systems will eventually syn- 
chronize. For IGR J 1 509 a recent polarimetric survey of new IPs 
has indeed revealed that the optical-red light is polarized at a few 
percent (~ 2%), suggesting a strong magn etic field (>10MG ) 
and a likely evolution towards synchronism (iPotter et al.ll2012h . 

4.2. Origin of pulses 

In most of the sampled sources the harmonics of the fundamental 
(up to 4to) are detected, indicating departures from strictly sinu- 
soidal pulse shapes. Non strictly sinusoidal pulses are seen also 
in the hard X-rays, where absorption is negligible. This suggests 
that the emitting regions depart from a symmetric dipole struc- 
ture. 
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V2069 Cyg is the one showing the most structured pulse in 
our sample (see figure|6), with a secondary maximum centered 
on a primary broad minimum seen at all energies. This indi- 
cates the presence of a secondary emitting pole. Furthermore, 
HRs in the 0.3-1 keV and 1-3 keV bands display a softening 
at spin minimum, indicating an additional contribution at this 
phase. V2069 Cyg also shows a BB component at soft energies 
(see Sect. 13.2. Il l whose contribution increases at the minimum 
of pulse (see table|6]l. Hence, a reprocessed component at the 
WD surface is present and is anti-phased with the optically thin 
emission from the PSR. 

On the other hand, 1GRJ1817 pulsates at a period that is 
half that found at optical wavelengths, indicating two symmet- 
ric poles contributing almost equally in the X-rays. The optical 
pulse is broad and overlaps the two X-ray maxima (figure|6) and 
hence likely originates from wide upper accretion regions above 
the two poles. 

For all sources the pulses are energy dependent, with ampli- 
tude decreasing at high energies. The softest band (0.3-1 keV) 
has pulse fractions even up to 99%, seen in IGR J0839, decreas- 
ing to a few percent in the harder band 5-15keV. The pulses 
harden at spin minimum in most sources as confirmed from the 
spectral analysis (see Sect. l3.2.2l l. This is a typical behaviour ob- 
served in IP type CVs and consistent with the accretion curtain 
scenario proposed bv lRosen et al.l (Il988l) . where the modulation 
is mainly caused by spin-dependent photoelectric absorption in 
an arc-shaped accretion curtain. The absorption is expected to 
be larger when viewing the curtain along the magnetic field lines 
(spin minimum). At these phases the softest regions are absorbed 
resulting in a spectral hardening. 

The curtain is also expected to contribute in the UV and op- 
tical light and hence the pulses at these wavelengths should be 
in phase with the X-rays. This is the case of 1GRJ0839, and 
IGR J 1509, though the optical light is dominated by the beat in 
the first system. A peculiar case is XSS J0056 that reveals opti- 
cal pulses changing phase with time: they are phase aligned with 
the X-rays only during the first 2/5th of observation. This may 
suggest a non-stationary contribution from two poles in the op- 
tical Hght. Also IGR J 17 19 and V2069 Cyg display anti-phased 
X-ray and optical pulses that could be due either to a more lumi- 
nous secondary pole or to the heated WD polar cap. For V2069 
Cyg the second option could be favored because of the detection 
of the BB component in its spectrum. 

4.3. X-ray emission components and energy budget 

The analysis of the average 0.3-100 keV spectra gives evi- 
dence of multi-temperature emission components as well as of 
complex intrinsic absorption. A low temperature (~0.1-0.7keV) 
plasma traces the regions closer to the WD surface. The oxygen 
H- and He-hke Unes detected in the RGS spectra of IGR J 17 19 
and IGR J 1509 give further evidence for the presence of low- 
temperature plasma. Their intensity ratio (Mauche 2002) in- 
dicates a temperature of ~0.2 keV and ~0. 1 keV, respectively. 
Since no forbidden component is detected in the He-like triplet, 
these lines form in high density regime (n > 10'^ cm"-^). The 
presence of Ne IX and Ne X and their ratio further indicate a 
temperature at ~0.3keV in IGR J1719. 

The intermediate temperature component is required in the 
PN spectral fits to account for the K^^ lines of Fe XXV (6.7 ke V) 
and Fe XXVI (6.97 keV). These are strong in the spectrum of 
XSS J0056 thus providing a satisfactory fit with only two mekal 
components. 



On the other hand, the high temperature component should 
originate in regions close to the shock. The wide range of values 
(~30-60keV) found in our sample would also imply a range of 
WD masses. The emitting hot regions have typically larger (~3 
times) emission measures (EM, normalizations in table |4) than 
those of the intermediate component, being that of IGR J 1650 ~5 
times. As EM oc the larger EMhot might indicate a larger 

emitting volume and hence, on average, a larger size. 

Hence, the multi-temperature spectra indicate a temperature 
gradient in the PSR though we are unable to map it with the 
present data. 

Different ai-e the cases of RX J0636 and V2069 Cyg that dis- 
play a non-negligible soft BB component with kTsB ~ 70 - 80 
eV. These t wo sources add them selves to the small group of soft 
X-ray IPs dAnzolin et al.ll20d8h increasing the roster to 15 sys- 
tems. The ratio of bolometric fluxes LBB/Lth. are ~ 0.85 for 
V2069 Cyg and -0.13 for RX J0636, respectively Therefore, 
this component represents a substantial fraction of the energy 
budget in V2069 Cyg, but not in RX J0636. The normalization 
of this optically thick component increases at spin minimum in 
both sources. It can be identified with the X-ray heated WD po- 
lar cap, whose projected area is larger at this phase. We estimate 
an average spot area in V2069 Cyg, aeB = 4.5 x 10'"* Dg^jp^ and 
hence a fractional WD area / ~ 7.6 x 10"^ ^900 ' ^'^^ ^wd = 
0.82 (see table|5]and also its notes for the lower limit estimate of 
distances). For RXJ0636 we estimate a^s - 6-3xlO'-'Dj|^p^ and 

f ~ 9 X lO^^Djj.^^^. These values are small, but comparable to 
what found in other soft X-ray IPs with similar BB temperatures. 

An additional component, ubiquitous in magnetic CVs, is the 
Kq, fluorescent line at 6.4 keV. It is detected in all sources of our 
sample with large EWs (~ 100-200 eV), arising from reflection of 
cold material. The lack of velocity shifts and the EW variability 
at pulse maximum and minimum in IGRJ0839 and IGR J1509 
suggest a likely origin at the WD surface. Such feature should 
be accompanied by a reflection continuum but our data do not 
allow to constrain this component. 

A lower limit to the mass accretion rate can be derived as- 
suming that the accretion luminosity is totally emitted in the X- 
rays: Laccr = GMqMwdRwd ~ Lbb + Lth, where Lth is the 
bolometric luminosity of the X-ray optically thin components 
and Lbb is the reprocessed X-ray luminosity, that is only in- 
cluded for V2069 Cyg and RX J0636. The estimated values are 
reported in column 3 of table |5] The secular mass accretion rates 
expected from angular momentum loss (AML) due to magnetic 
braking at the adop ted orbital periods, would be in the range 
of 10"^ - 10"'' Mo (iMcDermott & Taamlll989h . for the sources 
above the gap. Hence, unless they are located at a distance 10 
times larger than the estimated lower limits (very unlikely pos- 
sibly except for IGR J1830 and IGR J1719), the X-ray emission 
can not be solely representative of accretion luminosity. A sub- 
stantial fraction should be reprocessed and emitted in low en- 
ergy bands, as supported by the detection of periodic signals at 
the spin and/or bea t periods at UV and /or optical wavelengths, 
as firstly raised by iMukai et al.l (|1994|) . The only exception is 
XSS J0056 that falls in the period gap, for which the lower limit 
to the rate is compatible with the values e xpected either from 
magnetic braking or gravitational radiation (IWamerill995b . 

The spectral analysis has also shown the presence of dense 
(column densities up to ~ 10^^ cm"^) material partially ab- 
sorbing the X-ray emission. The spin phase dependence of ab- 
sorption indicates that the spin pulsations are mainly due to 
photoelectric absorption from pre-shock material in the accre- 
tion curtain. Furthermore, we found indication for the pres- 
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ence of absorbing ionized material in one source (IGRJ0839), 
whose signature is an OVII absorption edge detected at 0.76 keV. 
IGRJ 0839 adds as the third IP showing this feature (Mukai et al. 
I2QQU Ide Martino et al.li2008h . A warm absorber is commonly 
recognized in Low Mass X-ray Binaries (LMXB) through the 
presence of absorption edges and lines that have been related to 
ionized atmosphere above the accretion disc, giving rise to dips 
in close to edge-on systems (Diaz Trigo et al. 2006). The spectra 
at pulse maximum and minimum of this source do not provide 
strong evidence of changing in the edge parameters, suggesting 
that this material could be located within the binary, but not in 
the pre-shock magnetically confined accretion flow. A possible 
similarity to the LMXB dippers should be further investigated 
with phase-resolved spectroscopy along the binary period that 
we are unable to perform with the present data. 

4.4. Are hard X-ray IPs different? 

The hard X-ray sources analyzed here are confirmed as IP type 
CVs. Their spin periods are well within the range of previously 
known members of this group as depicted in figure[TT] where 
our sample is shown together with the current confirmed IPs. 
Orbital periods have been estimated combining X-ray and opti- 
cal results. Most of them are found above the 2-3 hr period gap 
except for XSS J0056 whose high asynchronism (P(j/Pn=0.05) 
might suggest it will never synchronize. The inferred spin-to- 
orbit period ratios are $ 0.12 consistent with most IPs. Hence, 
there is no striking difference between hard X-ray discovered IPs 
and other known IPs. Nonetheless, the orbital period distribution 
is now being also populated by the new hard systems at Pn k 5 
h. 

To inspect whether these selected INTEGRAL and Swift IPs 
possess massive WD we derived their masses using the maxi- 
mum temperature obtained from fitting the broad-band spectra 
with the PSR model and taking it as an estimate of the shock 
temperature (see table [Sj. These are obtained using Tshock - 
3/8GMwDA'mH/kRwD and adopting the WD mass-radius rela- 
tion from Nauenberg ( 1972). The mean WD mass of our sample 
is 0.86±0.07 Mq- Consistent res ults were a lso found in a recent 
study of 17 IPs made by Yuas a et alJ (12010 ). though they sug- 
gest a likely bias toward high masses (due to a high energy sam- 
pling). The average WD mass we found is also fully consistent 
with a r ecent determinat ion of CV primary masses of 0.8+0.2 
made bv lZorotovic et al.l (1201 II) . These authors also find that the 
WDs in CVs are substantially more massive than those in pre- 
CVs (0.67±0.21Mo) and than single WDs (^0.6 Mq). While the 
latter cases imply that the WDs increase their mass during the 
CV stage, the fact that IPs have primaries with masses similar to 
those of the whole CV population, suggests that the mass is not 
the driving parameter Therefore, we conclude that the hard X- 
ray detected CVs, though possibly biased toward high mass, do 
not show remarkable difference from the whole CV population. 

With the increase in the INTEGRAL and Swift exposures 
with time, most of the undiscovered IPs within ~1.5 kpc will 
be likely detected. This is the case of the faint RXJ0636, 
which only recently appeared in the new Swift/BAT catalogue 
(ICusumano et al.l 12010) . Therefore, the true population of IP 
type CVs is expected to be eventually unveiled by sensitive hard 
X-ray surveys of next generation mission. Moreover, as the new 
discovered IPs are being found in the previously poorly popu- 
lated range of long orbital periods (P>5 h), this will allow us 
a better understanding of the evol ution these systems through 
comparison with recent models, e.g lKnigge et alJ (1201 ih . 



5. Conclusions 

The main goal of this work is the characterization of a sample 
of 9 new hard X-ray selected CVs, to unambiguously identify 
them as magnetic systems of the IP type. The main results can 
be summarized as follow: 

- All sources are strong X-ray pulsators with periods between 
470-1 820 s and thus are confirmed CVs of the IP type. 

- All but two sources are spin dominated systems in the X- 
rays. IGRJ1817 and XSSJ0056, remain ambiguous cases. 
IGRJ0839, IGRJ1719, IGRJ1650 and RXJ0636 are beat 
dominated at optical wavelengths. 

- IGRJ0839, V2069Cyg and IGRJ16500 are disc-fed accre- 
tors, while IGR J1830, IGR J1719, RX J0636 and IGR J1509 
display a disc-overflow accretion mode. IGRJ1817 and 
XSS J0056 could be either stream-fed or disc-fed systems. 

- We also estimated the orbital periods (except for IGR J1650), 
all found to be above the 2-3 h period gap. XSS J0056 is 
the only exception, lying in the gap. We derived first esti- 
mates forIGRJ0839 (8 ± 1 h), IGR J 1830 (4.2 + 0.2 h), and 
IGRJ1817 (6.3-8.7 h). 

- The amplitude of the X-ray pulse was always found to de- 
crease with energy due to complex absorption. 

- We identified two or three optically thin components with 
characteristic temperatures in the ranges: 0.1-0.7, 4-10 and 
30-60 keV (signature of a temperature gradient in the post- 
shock region). Three are present in six sources. 

- V2069 Cyg and RX J0636 are instead found to possess a soft 
X-ray optically thick component at kT^^ ~80 eV. This in- 
creases to 15 the current roster of soft X-ray IPs, confirming 
a relatively large (~ 30%) incidence. 

- A strong K(j, Fe line at 6.4 keV is present in all sources likely 
originating at the WD surface from the reflection of X-rays 
from cold material. 

- The spectrum of IGRJ0839 reveals an absorption edge at 
0.76 keV from OVII. This increases to three the IPs were a 
warm absorber is found. 

- The WD masses of our sample are found between 0.74 and 
0.96 Mq, with an average mass 0.86 + 0.07 Mq, suggesting 
no striking difference with the whole CV population. 
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